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ABSTRACT 

Hierarchical linear models distinguish between the 
individual and the group levels of data. Hence, they are often 
referred to as multilevel models. It is easiest to think of 
hierarchical linear models as special regressio: models that allow 
simultaneous investigation of the respective roles that individual 
and group characteristics play in the attainment of treatment goals. 
A hierarchical linear model can be fit to the data to investigate the 
contribution of individual characteristics to young people's end of 
treatment attainment and to determine how these relations are 
influenced by group characteristics • To simplify, people think of 
data as consisting of only two levels: individuals nested within 
treatment groups; and groups. The data are thus described by two 
models: the first model is based on individual data. It describes the 
prediction of individual outcome from individual characteristics and 
is referred to as the within-unity, or unit level mode. The second 
model expresses the variability in regression coefficients as a 
function of group-level variables. This model is referred to as the 
group-level, or between-unic mode. The majority of research studies 
on group therapy have ignored the group effect. Given the fact group 
characteristics are the prime reason for selecting this treatment 
modal ity, it appears important to include these characteristics in 
any model for treatment effectiveness. (ABL) 
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An Application of Hierarchical Linear Modeling 
to Group Research 



The purpose of this discussion is to introduce the 
group researcher to statistical methods better suited to 
the study of group process than traditional analysis of 
variance and multiple regression. These methods are known 
by the name hierarchical linear models. Hierarchical 
linear models distinguish between the individual and the 
group levels of data. Hence, they are often referred to 
as multilevel models. It is easiest to think of 
hierarchical linear models as special regression models 
that allow simultaneous investigation of the respective 
roles that individual and group characteristics play in 
the attainment of treatment goals. The particular 
individual characteristics of interest will depend on the 
nature of the treatment group. Group characteristics of 
interest typically pertain to leader qualities and 
contextual factors that describe the effects of the group 
members on one another. 

Despite the nearly universal agreement that individual 
and group characteristics are important influences on 
individual outcomes, the literature coes not reflect much 
empirical interest in individual characteristics and 
ignores group characte*-istics. The lack of empirical 
interest in group characteristics is at least partially 
due to the difficulty in using traditional statistical 
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methods to study hypotheses that are implicitly 
multilevel — i.e., when key independent variables, such as 
group characteristics, are typically measured at a higher 
level of aggregation than the outcome variable of interest 
(individual attainment on some attribute of interest). 

Common individual characteristics such as gender and 
ethmcity have occasionally been included as factors in 
group outcome studies. These studies have used factorial 
analysis of variance, which permits the investigation of 
differential treatment effectiveness for classes of 
individuals with certain characteristics. By far the more 
frequently used procedure has been one-way analysis of 
variance, which simply permits the investigation of 
systematic c.fferences among treatment means. Individual 
characteristics are treated as a nuisance, or source of 
error in the data. Studies using traditional analysis of 
variance techniques are typically characterized as a 
comparison between a treatment group and a control group, 
or between one therapeutic orientation versus an 
alternative one. Both designs ignore the fact that 
individuals may differ in their response to similar 
treatments as a consequence of the unique developmental 
ecology of each group in which the treatment is 
implemented. Traditionally, such individual and 
group-level variation is seen as a nuisance, in fact, the 
traditional definition of error variance implies that 
indiviauals within a group provide independent responses, 
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This assumption is usually untenable in group research 
where treatments are administered to groups of subjects. 
As a result, the standard error of the treatment effect 
will typically be underestimated, leading to 
overestimation of the statistical significance of 
treatment effects. Alternatively, the group instead of 
the individual can serve as a unit of analysis. 
Unfortunately, this strategy lacks power for detecting 
treatment differences. Moreover, it prevents differential 
treatment effectiveness for individuals with different 
characteristics. It also prevents the use of individual 
characteristics, such as pretest scores, as covariates. 
Since covariates can, in principle, dramatically enhance 
the statistical power of such a study (Porter & 
Raudenbush, 1987), this shortcoming is a serious on^ 
(Raudenbush & Bryk, 1988). In both the individual and 
group level approaches, the multilevel structure of the 
data is lost, i.e., neither analysis recognizes that 
individuals are nested within groups. Not surprisingly, 
the opportunity is also lost to investigate the dynamic 
interplay of individual and group characteristics in the 
attainment of treatment goals. 

Failure to investigate the interaction of individual 
and group characteristics in the attainment of treatment 
goals has been cited as a serious shortcoming in an 
earlier review of group outcome studies (Bedner & Kaul , 
1978). In general, this review suggests that although 

b 
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research justifies group treatment, and provides a 
stimulus for further research, one must be aware of what 
is not being studied. Specifically many of the primary 
and unique variables of group treatment are not being 
subjected to empirical test. 

The inability to adequately test the primary and 
unique variables of group treatment with traditional 
linear models has not prevented the conceptualization of 
these variables. Three noteworthy sources of information 
on the mechanisms of groups are Corsini and Rosenberg 
(1955), Hill (1957), and Yalom (1975). These works 
identify similar variables that may be subsumed within 
three higher-order concepts. First, group members may 
improve as a consequence of learning based on their 
participation in, and evaluation of, a developing social 
microcosm. Second, psychological growth may result from 
soci al 1 earni ng processes based on i nterpersonal f eedbacK 
and consensual validation. And third, individuals may 
profit from the reciprocal opportunities to be both 
helpers and helpees in group settings (Bedner & Kaul , 
1978). Without an adeauate methodology to test the 
validity of these concepts, the likelihood of any 
integrated understanding of the uniqueness o"^ grouo 
treatments seem low. 

The development of adequate computer software based on 
suitable mathematical models for data with a multilevel 
structure occurred in the mid 1980's. The reasons for 




this development were concern about the inferences 
traditional methodologies yielded when data occur on more 
than one level, and the lack of compatibility in the 
conceptualization of multilevel data and the statistical 
analysis of these data. 

Educational researchers (Aitkin & Longford, 1986; 
Goldstein, 1986,1987; Raudenbush & Bryk, 1986; Willms, 
1987) who focus on the growth of student knowledge and 
skill occurring within educational organizations have been 
most active in analyses of multilevel data. The 
discussion and application of hierarchical linear models 
which follow are based on adapting innovations from 
educational research to group research. 

Consider the following multi-level data from a group 
counseling situation. Therapy groups are formed at 
several juvenile probation departments for tne purpose of 
enhancing juvenile offenders' level of social development. 
An individual measure of social development is obtained at 
the outset and termination of therapy. Offenders are 
assigned to weekly group counseling as a condition of 
their probationary status for a period of one year. The 
offenders dif'Fer in age and severity of the offenses 
committed. Measurements are made on group characteristics 
such as experience of the group leader, and the average 
interpersonal feedback provided from group members. 

A hierarchical linear model can be fit to the da^ca to 
investigate the contribution of individual characteristics 
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to young people's end of treatment attainment and to 
determine how these relations are influenced by group 
characteristics. To simplify, we think of the data as 
consisting of only two levels: individuals nested within 
treatment groups; and groups. The data are thus described 
by two models: the first model is based on individual 
data. It describes the prediction of individual outcome 
from individual characteristics and is referred to as the 
within-unit, or unit level model. In the context of the 
present example, this model states that post-test social 
development is predicted from pre-test social development, 
age of offender, and the severity of the committed 
offense. If one were simply studying individuals within a 
single group, standard regression techniques could be used 
to estimate the coefficients of this first equation. 
However, considering the multilevel structure of the data, 
it is reasonable to presume that the regression 
coefficients of the first equation vary across trea^.^ent 
groups. In fact, it is this variation which is of 
particular interest to group researchers. If groups play 
a role in indivi::'u:. outcome, then group characteristics 
will predict variation in the coefficients of model 1 
across groups. As was stated before, the data from this 
example are described by two models. In fact, the second 
model expresses the variability in regression coefficients 
as a function of groups-level variables. This model is 
referred to as the group-level, or between-umt model. In 
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this model, the parameters (intercepts and slopes) derived 
from the within-unit model for each group are expressed as 
a function of group characteristics. In the context of 
the present example, these characteristics are the 
experience of the group leader and the average 
interpersonal feedback provided by group members. 

It is perhaps easiest to think of the two models as 
two separate stages in an analysis strategy. At the first 
stage, a regression analysis is performed in each 
treatment group separately. The coefficients from these 
analyses are saved and then become the dependent measures 
in the second stage of the analysis. In this latter 
stage, the estimated coefficients from the first stage are 
regressed on the group characteristics of leader 
experience and average interpersonal feedback from group 
members. While it is convenient to think of the analysis 
as consisting of two separate stages, in truth both models 
are fit simultaneously. To formalize our discussion 
somewhat, it is necessary to introduce some notation. The 
within-unit, individual level equation is: 

Yij = Bjo -h Sj/Xyji -h Bj2Xijz + BjaXijs + Rij 



where 

Yjj is the post-treatment social development score for 
individual 7 (7 =1 , . .nj ) in group j ( j = 1 . . .k) • 
Bjo is the constant intercept for group j. Bjo 

erJc i) 
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describes how well an individual with average 
individual pretest, age, and severity of offense 
performs on the social development posttest measure. 
In other words, Bjo is an estimate of the expected score 
in group j after controlling for variability in 
pretest score, age, and severity of offense at the 
level of the individual. 

represents the pretest social development score 
for individual 7 in group j. 

Bj 1 is a regression coefficient that describes the 
effect of the pretest on the outcome measure of social 
development in group j, holding constant the effect of age 
and severity of offense, 

Xij2 represents the age of individual 7 in group j, 
Bj 2 is a regression coefficient for the effect of age 
on the outcome measure of social development in group j 
holding constant the effect of pretest and severity of 
offense , 

Xij3 represents the severity of offense of individual 
7 in group j\ 

Bj3 is a regression coefficient for the effect of 
severity of offense on the outcome measure of social 
development in group j holding constant the effect of 
pretest and age. 

Ri j is the unique contribution of irdividual i in 
group J, Rfj is the extent to which the outcome for 
subject 7 in group j is not fully explained by preceding 

10 
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terms in the model. 

From the previous description, we can see that it is 
the BjS {Bjo, Bj 1 , Bj 2 , Bjs) that are presumed to vary 
across groups, and that are to be predicted by the 
group-level characteristics. Because these parameters 
vary, they are known as random coefficients, or random 
effects. 

The between-unit model is a system of equations based 
on the group level data; it is a model for the random 
coefficients from the within-unit model and is expressed 
as : 

Bjo = Ooo + O01D1J + Oo2D2j + Ujo 
Bjp = Opo + OpiDij + 0p2D2j + Ujp 

where 

Bjp is the pth coefficient f'^om the unit level 
equation for group j. To make this more concrete, the 
between unit model will be described specifically for 
Bjo, the intercept of the within-unit equation. 

Bj 0 is the intercept for group j from the within-unit 
equation. As such, it describes the expected post- test 
for an individual from group j with average pretest, age, 
and severity of offense. 

Ooo is the constant intercept from the group level 
model. Thus, Ooo is the expected value of the unit level 
intercept (Ej (Bjo) when experience and interpersonal 

11 
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feedback are held constant. Thus, Ooo represents the 
grand mean across groups in their adjusted level of 
performance. 

Dij is the score for leader experience in group J. 

001 is a regression coefficient describing the effect 
of leader experience on expected group attainment, holding 
constant the effect of average group interpersonal 
feedback to others. 

D2J is the average interpersonal feedback score for 
group J. 

do2 is a regression coefficient describing the effect 
of average group interpersonal feedback on expected group 
attainment, holding constant leader experience. 

Ujo is the unique contribution of group, i.e. Uj o 
describes the extent to which Bjo is not exactly 
reproduced by the group level model. 

As mentioned, additional group-level regression 
equations are included in the between unit model. One 
equation is provided for each coefficient from the 
within-^unit equation. This provides a multivariate 
formulation for examining the effects of between-group 
variables (e.g., leader characteristics) on within-group 
relations (e.g., the age-attainment relation). 

This example was based on a data structure with two 
levels. It is possible to fit data with higher level 
structures. For example, it is possible that there are 
relevant differences across different juvenile probation 
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departments. Department could be incorporated as a third 
level in the model. In this case the structure of the 
data would be individuals within groups, groups within 
departments, and departments. Hierarchical linear 
mode ling can al so be used i n group process stud i es . Here 
the data would be structured as time within individuals, 
individuals within groups, and groups. Once one has a 
multilevel mindset in regard to date., it is easy to see 
the complex structure of data that was previously 
conceptualized as single level. 

The computational software for multilevel models is 
available from several groups of investigators (Bryk, 
Raudenbush, Seltzer, & Congdon, 1986; DeLeeuw & Kreft, 
1986; Goldstein, 1986; Longford, 1987; Mason, Wong, & 
Entwisle, 1984). The statistical theories upon which 
these snproaches are based and therefore the statistical 
properties of their results are, for most practical 
purposes, identical. All use iterative procedures to 
compute maximum likelihood estimates of variances and 
covariances. These estimates are computed by different 
numerical approaches which need not concern us here. 
Interested readers can see Raudenbush and Bryk (1988). 
The estimation procedures share five basic properties: 1) 
the effects of the sampling variance are minimized, 2) the 
covariation among the B coefficients are taken into 
account, 3) total variance in parameter estimates can be 
partitioned into true parameter variance and sampling 

13 
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variance, 4) parameter covariation among the Ss is 
provided, and 5) the estimators of Bjk are corrected for 
unrel i abi 1 i ty . 

As can be seen from this example and discussion, the 
implementation of hierarchical linear modeling greatly 
expands the range of methods for investigating groups. It 
also has the potential to extend conceptualization about 
relevant individual and group characteristics. It must be 
pointed out, however, that the design and data 
requirements to produce an adequate specification of the 
within-unit and between-unit equations are demanding. 
Identification of relevant individual and group 
characteristics is the first design requirement. A second 
requirement concerns obtaining a valid and reliable 
measure of attainment of the treatment goals at the 
individual level. Third, a pre-test measure should be 
built into the analysis to ensure that effects which any 
model may attribute to higher-level units, such as groups, 
are not simply a consequence of pre-existing differences 
among individuals within groups (Garner, 1989; Hauser, 
1970). A fourth design point concerns sufficiently large 
samples at the group level to provide accurate estimates 
of group level parameters (Goldstein, 1987). Ideally, 
about 50 groups would provide that variation. 

In summary, a methodology has been presented for the 
adequate assessment of the individual attainment of 
treatment goals that occur within groups. The majority of 
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previous research studies on group therapy have ignored 
the group effect. Given the fact that group 
characteristics are the prime reason for selecting this 
treatment modality, it appears important to include these 
characteristics in any model for treatment effectiveness. 
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